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The CLEO Collaboration recently reported the observation of the B" — » D*~'pn mode at a rate only 
a factor of 4-5 lower than S" — > D*~ p^ and D'^-k^ . We try to understand this with a factorization 
approach of current produced nucleon pairs. The baryon weak vector current form factors are 
related by isospin rotation to nucleon electromagnetic form factors. By using G^/* measured from 
e'^e~ — * NN and pp — > e"'"e~ processes, assuming factorization of the B^ -^ D*~ transition and 
the pn pair production, we are able to account for up to 60% of the observed rate. The remainder 
is argued to arise from the axial current. The model is then applied to B decays to other mesons 
plus pn modes and D*~ plus strange baryon pairs. 



PACS numbers; 13.25.Hw, 13.40. Gp, 14.20.Dh 
I. INTRODUCTION 

Following a suggestion by Dunietz |lj that B -^ 
DNNX decays could be sizable, the CLEO Collabora- 
tion has recently reported the first observation of such 
modes |2J: 



Br(B° ^ D*-pn) 
Br(B° ^ D*-ppTT+) 



(14.5^3:0 ± 2.7) X 10-^ (1) 
(6.5ti:2±l-0)xl0-^ (2) 



Although the decay final states are three or four-body, 
they are only a few times below corresponding two-body 
mesonic modes H such as 



Br(S" 
Br(S" 



D* 



= (6.8 ±3.4) X 10 



-3 



D*-7r+) = (2.76 ±0.21) x 10" 



(3) 
(4) 



Since D*~ creation carries away much energy, the ob- 
served large rate of i?° -^ D*~pn supports the suggestion 
that enhanced baryon production is favored by reduced 
energy release on the baryon side B. Thus, given the 
large rate oi B -^ rj' + Xs decay where p„' > 2 GeV [^ , 
the B -^ T]'Ap decay may be sizable [Q compared to 
charmless two body baryonic modes. Similar argument 
holds for B -^ ^Kp as implied by i? -^ 7 + Xg. Since 
A — > pi: decay automatically provides spin information, 
the observation of such charmless three (or more) body 
baryonic modes involving A baryons allows for a search 
program for triple-product type of CP and T violating 
effects. With this in mind, a better understanding of the 
B^ — > D*~pn decay is not only worthwhile in its own 
right, it can also serve as an important first step towards 
a more ambitious project on charmless baryonic modes. 
In order to account for the proximity of rates shown 
in Eqs. (|^) and (^, B^ — > D*^pn was seen through 
a simple pole model M as occurring in two steps of an 
underlying b -^ cmJ transition, i.e. B -^ D*~h followed 
hy h ^ NN, where h stands for the off-shell 77+ and 
/9+ mesons, as illustrated in Fig. ||(a). Taking this as a 
Feynman diagram, the decay amplitude of i?" — *■ D*~pn 
can be written as 
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where the B^ -^ D*^ transition matrix element is 

2V{q^) 



v{pn), (5) 
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a 
AuapPEPD' 



ms + niD* 



9i^u 



<itj.qu 



[niB +mD')Ai{q^) 



+i ipB+PD')f, 



-% 



A2{q^ 



ms + riiD 



-Qu 



-i 2to£). 



'-Ao(q' 



= D*i 



(6) 



q = Pb — Pd* = Pp + Pn is the momentum transfer (so 
t — q^ is nothing but the pn pair mass), eo* is the po- 
larization of the D*~ meson, and gf and ^2 ^^^ 
respectively the vector (Dirac) and tensor (Pauli) cou- 
pling constants of the p-meson to the nucleons. 

In Eq. (H), factorization of the first vertex is justi- 
fied to some extent, and is expressed through the effec- 
tive coefficient ai. For example, in "naive factorization", 
ai = ci + C2/Nc where Nc is the number of quark col- 
ors and ci,2 are Wilson coefficients, but in general the 
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FIG. 1. (a) The simple pole model and (b) the factorized 
current model for describing B° —* D*~pn. 



effective coefBcient ai can be extracted from experimen- 
tal data on i3° -^ D*~ p'^ decay. The problem of the 
pole model picture is with single p dominance and the 
q^ dependence of the strong interaction pNN vertex. 
That is, there is no reason why the ud weak current 
only generates a p meson, which then propagates and 
generates the baryon pair, as seen from the second and 
third lines of Eq. (^) . Quantitative results based on this 
model is highly unreliable, which is further aggravated 
by our ignorance of time-like meson- nucleon couplings. 
Since our knowledge of the coupling gf are based on 
low Q^{= —q^) space- like processes, they cannot be ex- 
pected to give reliable quantitative results for time-like 
processes at higher energies. 

In this paper we turn to a different approach by propos- 
ing a generalized factorization of current produced pn 
pairs. The three-body decay is seen as generated by two 
factorized weak currents (linked by a VF-boson), where 
one current converts B'^ into D*~ and the other creates 
the pn pair, as shown in Fig. ^(b). In this way, having 
factorized the 5° -^ D*^ transition, we concentrate on 
the weak current production of baryon pairs. 

It is well known that the vector portion of the weak 
current and the isovector component of the electromag- 
netic (em) current form an isotriplet. Thus, information 
on the nucleon em form factors, for which much data ex- 
ist in both the space-like region |q] (e.g. ep — > ep), and, 
of particular interest to us, the time- like region |7[-[lC|] 
(e.g. e^e~ — *■ NN and pp —^ e~^e~), can be transferred 
to the nucleon weak form factors by a simple isospin ro- 
tation. The total decay amplitude that comes from the 
vector portion of the weak current can be written down 
unambiguously and the portion of the branching fraction 
that comes solely from the weak vector current can be 
readily obtained once the form factors are given. We 
find that the vector current can account for 50%-60% of 
the observed rate of Eq. (1^). Since this analysis involves 
just the factorization hypothesis but is otherwise based 
on data, it is rather robust. 

Although the axial vector and vector current contribu- 
tions interfere in the decay amplitude, the interference 
vanishes when one sums over spin and integrates over 
phase space. The total rate is therefore a simple sum of 
the contributions from the vector and axial vector por- 
tions of the weak nucleon form factor. Like the vector 
case, we could in principle obtain the axial vector con- 
tribution if the nucleon form factors of the axial current 
were known. Unfortunately, the time-like data is still 
lacking, hence the contribution from this part remains 
undetermined. In spite of this, a rough estimate can still 
be given, which seems to be the right amount. We point 
out, however, that information on the time- like nucleon 
axial form factor could in fact be obtained in the future 
via the B° — » D*~pn decay data. One just has to sepa- 
rate the axial vector contribution from the vector part. 

This paper is organized as follows: in the next section 
we lay out the factorization assumption that allows us 



to relate the vector current contribution to the nucleon 
em form factors, where one enjoys abundance of exper- 
imental data. We are then able to compute the vector 
current contribution to B'^ —> D*~pn. The axial vec- 
tor contribution is also estimated. 
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In Sec. 

trate the power of our data based approach, we briefly 
introduce an improved pole model approach and stress 
the need for improved measurements of neutron em form 
factors. Finally, we apply our analysis to B^ -^ D^*^^pn 
and i?" — > D~pn and other baryonic modes containing 
strangeness, and conclude in the last section. 



II. FACTORIZATION APPROACH AND 
NUCLEON FORM FACTOR DATA 



As illustrated in Fig. pj(b), we generalize factorization 
from two-body to three-body decay processes, 

{D*-pn\neff\B") = ^VudV:,a,{D*-\V^^ - A^\B") 



x{pn\Vf,-A^\Q). 



(7) 



The first matrix element is as before, but for the sec- 
ond, the nucleon pair is viewed as directly created by the 
current. The vector part can be expressed as 



F^ (t) 
{pn\V^\0) = u{pp) \ F^{t)j^ + i^fAi^^^q- k „(p,_^) 



2niN 



(8) 



where ttijv is the nucleon mass, t = {pp + Pn)"^ = q^, and 
F]*^ are the weak nucleon form factors. Likewise, the 
weak axial current A^ = A}^ + iA^ matrix element is 

{pn\A^\Q) = u (pp) \ gA (t) 7^ + -^ g^ f ^^ " (^«) ' 

L 2 TUN J 

(9) 

where gA{t) is the axial form factor, and hA(t) is referred 
to as the induced pseudoscalar form factor. Expressions 
for the space- like processes are similar. The weak nucleon 
form factor F^ (t) and the axial form factor gA{t) are 
normalized at t = H 

Ff'(O) = 1, 5a(0) =gA = 1-2670 ±0.0035. (10) 



A. Isospin Relation and Nucleon Form Factors 



It is well known that the photon field A^ contains W^, 
which forms a weak isotriplet with W^'^. Therefore the 
currents they couple to also form an isotriplet, and are 
related by an isospin transformation. For the nucleon 
system, the strong isospin symmetry coincides with the 
weak isospin symmetry of the weak and em currents. The 



weak vector form factors are therefore related to isovector 
electromagnetic (em) forra factors. 

The matrix element {N{j)')N{p)\Jli"^\G) for the em 
current can be expressed as 



{N{p')N{p)\j;^^\Q) 



(11) 



Similar expressions can be obtained for space-like pro- 
cesses. The quantities Fi{t), F2{t) are respectively the 
Dirac and Pauli form factors, normalized at i = as 

i^f(0) = l, F«(0)=0, F|(0) = ..p, F2"(0) = «;„, 

(12) 

with Kp („) the proton (neutron) anomalous magnetic 
moment in nuclear magneton units. The experimental 
data is usually given in terms of the Sachs form factors, 
which are related to Fi and F2 through 
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(13) 



One clearly has Gm = Ge at threshold t = 4,171%, while 
at i = we have 

G^(0) = 1, G^(0) = 0, Gl,{0) - 1 + «„ GIj{0) = Kr,. 

(14) 

The isospin decomposition of the em current is given 
by the following definitions 
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{F[±Fn, (* = 1,2), 



(15) 



where Ff and F" are the isoscalar and isovector decom- 
positions of the form factors, respectively. The fact that 
the isovector component of the em current, together with 
the vector portion of the charged weak currents, form an 
isotriplet is manifested by 



2F^it) = F.rit), 



z = l,2, 



(16) 



the factor 2 coming from the definition of -^{'2^ (i) in 
Eq. (|l|). For example, from Eqs. (0), (|l|) and ^ we 
have2Fl'{0)=F^{0). 

With Eq. (|6|) and the Gordon decomposition, we can 
put the three-body S" — ^ D*~pn decay amplitude of 
Eq. (0) into the following form 



iMv = -i^VudV,laie*j^. 
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(17) 



where V indicates that the nucleon pair is generated 
by the vector current. The terms proportional to q^, in 
Eq. (pi) vanish in the limit of equal proton and neutron 
mass. For completeness, the amplitude for nucleons gen- 
erated by the axial vector current is given by 
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V [Pn] , (18) 



where (£'*"|V^ - ^mI^°) is given in Eq. (|). 

The two amplitudes Ma and My interfere, since 

Y,'2Re{MAM*v)^32G%\Vudf\V,bfal 

xgA{t)Vit)A,{t)(^Glj{t)-Glj{t)) 

X {PB ■Pp){PD* -Pn) - {PD* ■Pp){PB ■ Pn) (19) 

is in general non- vanishing. The summation is per- 
formed over the nucleon spins and D*~ polarization. The 
three-body phase space is described by the two inde- 
pendent variables rripjj = t = {pp + pnY a-^d rri]-!,^ = 
{pd* +Pn)^- The interference term is antisymmetric with 
respect to exchange of Pp and Pn- For given t, as we in- 
tegrate over the kinematically allowed rn?^,,^, each value 
of ^ 2Re{MAMy) from a given pair of Pp, pn would be 
cancelled by those from the exchanged pair. The inter- 
ference term thus contributes nothing to the total three 
body decay rate F, and the final result is a simple sum 
of the contribution from My and M.a separately. 

It is interesting to note that, although the effect of 
^ 2i?e(A^A-A^y) does not show up in the decay rate, 
we can nevertheless construct an asymmetry ratio mea- 
surable based on the antisymmetric nature of this quan- 
tity, to extract the time- like information of gA{i) from 
B° —f D*^pn data. The information of V{t), Ai{t) and 
G^{t) in Eq. ( p^ ) can be found by other means, and 
the overall factors ai and \Vud\, \Vcb\ would cancel in the 
asymmetry ratio. 



B. Form Factor Data and Perturbative QCD 

Much data has been accumulated for the nucleon em 
form factors, which turns the vector portion of the de- 
cay amplitude iAiy expressed in Eq. ( |l7|) into something 
that we can handle. The branching fraction can be read- 
ily obtained once the nucleon em form factors are given, 
thanks to the isospin relation of weak vector and em 
currents. It is important, however, to make sure that 
the form factors satisfy the constraint from perturbative 
QCD (PQCD). 

The so-called quark counting rules ||ll| give the lead- 
ing power in the large-|t| fall-off of the form factor Fi{t) 
by counting the number of gluon exchanges which are 
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FIG. 2. Time-like proton magnetic form factor data, fitted 
(solid) by Eq. (E2|) with the parameters given in Eq. (|24|). The 
other (dash) curve is discussed in Sec. III. 



necessary to distribute the large photon momentum to 
all constituents. Since helicity-flip leads to an extra \/t 
factor in F2(t), one finds, in the limit |i| —> oo 



F.{t) ^ (|i|)-(^+i) 



m4 



i^ 1,2 



(20) 



where /3 is the /3-function of QCD to one loop, and Qq ~ 
Aqcd — 0.3 GeV. We note that 7 depends weakly on the 
number of flavors; for three flavors 7 = 2.148. 

The asymptotic form given in Eq. ( |20| ) has been con- 
firmed by many experimental measurements of Gm — 
Fi + F2 over a wide range of momentum transfers in 
the space-like region. The asymptotic behavior for G^.j 
also seems to hold in the time-like region, as reported 
by the Fermilab E760 experiment ||] for 8.9 GeV^ < 
t < 13 GeV^. Another Fermilab experiment, E835, 
has recently reported G\.j for momentum transfer up to 
~ 14.4 GeV^, and gives the empirical fit ||]: 
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which agrees well with the asymptotic form in Eq. (pof). 
By exploiting the relation in Eq. (nSh , the combination 
2{Fl^F}i) in Eq. ^ can be replaced by G\^-G\i which 
is composed of measurable quantities. Similar replace- 
ment can also be made for _F2, which is a combination of 
G\j^—G\ and G^ — G^. Most time-like data for the mag- 
netic form factors, however, are extracted by assuming 
either \Ge\ ~ |Gj\/| or \Ge\ = in the region of momen- 
tum transfer explored. Since Gm — Ge = (1 — t/Amj^)F2 
clearly vanishes at threshold, the absence of clear devia- 
tion from this assumption in extracting data implies the 
contribution of F2 is negligible even for t far beyond the 
threshold, which is consistent with the prediction from 
QCD. In fact, by assuming \Ge\ — \Gm\ in extracting 



FIG. 3. Time-like neutron magnetic form factor, where 
solid (long dash) line is fitted by Eq. (p3) with parame- 
ters given in Eq. (051) (Eq. (p6|)), for data extracted with 
\Ge\ = |Ga/| HGeI = 0) assumption. The third (dash) curve 
is discussed in Sec. III. 



G]\i from data, the information on F2 is lost. In our cal- 
culation we concentrate on the part of Eq. (HTh which 



contains F^ 



Fl 



The contribution from Fg can be de- 



termined only when Gm and Ge can be separated from 
data with better angular resolution. 

We take \Gm\ in the following form to make a phe- 
nomenological fit of the experimental data |0-UQ| : 



|GS,(0l = (|iH-f + 5 + ^ + 5 
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where the leading power and log factor are from Eq. (& 
and the fewer parameters for G^lj reflects the scarcer 
amount of neutron data. We find the best fit values 



2^3 

and 



a;i= 420.96 GeV* 
X2= -10485.50 GeV 
106390.97 GeV^ 



Xi 



-433916.61 GeV' 
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X5 = 613780.15 GeV' 



(24) 



j/i = 236.69 GeV* 



j/2 = -579.51 GeV^ (25) 



where the x^ per degree of freedom of the fits are 1.47, 
0.41 for \G^m\^ I'-'mI' respectively. We show in Figs. § 
and H the fitted data together with the best fit curves 
given by Eq. ( p^ ) and (|2^) with the parameters above. 

It was pointed out in Ref. [g| that the data supports 
|G^| = as well. We therefore perform the fit for the 
neutron magnetic form factor to the data that is ex- 
tracted under the assumption of |G^| = 0. Since the 



number of data points is small, a two-parameter fit as in 
Eq. (E3) would still suffice. The best fit values are 



yi = 292.62 GeV" 



2/2 = -735.73 GeV^ , (26) 



giving x^/d.o.f. = 0.39, which is a little smaller than 
the previous fit, and the fit is plotted as the long-dashed 
line in Fig. ^. To the eye, the data might slightly prefer 
the second fit, especially the t — 6 GeV^ data point. 
However, it should be clear that more data is needed to 
distinguish between these two cases. 

We note that there is a sign difference between G^ 
and G%[ in the space-like region, since from Eq. (Ilj) 

GIM ■ ^ ■ " "■ ■ 



1 



.p > and GIj{0) 



< 0. Analytic- 
ity implies continuity at infinity between space-like and 
time-like t P4 I regions. Hence time-like magnetic form 
factors are expected to behave like space-like magnetic 
form factors, i.e. real and positive for the proton, but 
negative for the neutron. 

For large t, QCD predicts the magnetic form factors 
to be real |ll| , with the neutron form factor weaker than 
the proton case |lq|. According to QCD sum rules p4[ , 
asymptotically one expects G^,f/G^j^ ~ Qd/Qu = —0.5. 
We can readily check our fits against these: for Gm 
fitted to data extracted by assuming \Ge\ — \Gm\ for 
both neutron and proton magnetic form factors, we have 

GIj/GIj = -yi/xi = -236.69/420.96 == -0.56 0.5,. 

For Gm fitted to the proton data extracted with the as- 
sumption \Ge\ — \Gm\ but the neutron data extracted 
assuming \Ge\ — 0, we have G'^/G\j = —yi/xi = 
-292.62/420.96 = -0.70, slightly larger than -0.5. Nu- 
cleon form factors have also been analyzed from nega- 
tive to positive t with dispersion relations. The phase 
of the proton magnetic form factor turns out to be 
~ 27r, hence the proton magnetic form factor is real and 
positive as expected asymptotically, but already from 
t > 4 GeV^ (l5|jl|] onwards. 
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for the BSW model. The subscript V serves as a re- 
minder that this is the result from the vector portion 
of the weak current alone. The upper and lower errors 
correspond respectively to the maximum and minimum 
of the branching fraction evaluated by scanning through 
both x^ < xLn + 1 Of \Gl,\ and \Gl,\ fits. 

For proton data extracted by assuming \G^^\ 
but neutron data assuming |G^| — 0, we have 
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(29) 
(30) 



The larger values for this second set of branching frac- 
tions can be understood qualitatively from Fig. pi where 
the curve fitted to data assuming \Ge\ = is higher than 
the one fitted assuming \Ge\ — \Gm\- Since the neutron 
magnetic form factor is negative in the time-like region, 
the quantity G\j — G\.j is larger if G^j gets more negative, 
and the branching fraction becomes larger. 

Comparing with the central value of the measured 
Br(SO -^ D*-pn) = (14.5+^-o ± 2.7) x 10"*, our 
LF (BSW) model results contribute 50 (60)% for the first 
set and 62 (75)% for the second. If the naive factorization 
vaue for oi is used, the results are close to the experimen- 
tal central value, that is Br^'"^^^^^ = 12.51 (13.83) x 
10-4 for the first set and 15.69 (17.36) x 10"'' for the 
second set. We plot in Fig. the vector current induced 
differential decay rate dTy [B^ -^ D*~ pnj / dq^ , for both 
the LF and the BSW models. The lower two curves are 
from the approach of next section. As seen also from 
Eqs. ([27|-|30|), the LF form factor model gives results that 



C. Results for B° -^ D*' 



pn 



Before using data and the nucleon form factor rela- 
tions to compute the decay rate, we need to specify the 
value of ai to be used. Since Eq. (O) depends on the 
product of fli and B -^ D* form factors ||l^, we take 
a phenomenological approach and use the value of oi 
extracted from the two-body mode B -^ D*^ p^ , i.e. 
qBSW = o.86± 0.21 ±0.07 and o^f = 0.74±0.18±0.06for 
Bauer-Stech-Wirbel (BSW) |l|,|l| and light-front (LF) 
form factor models [Q. The CKM matrix elements Vud 
and Vcb are taken to be 0.9757, 0.039 respectively. 

For both proton and neutron data extracted by assum- 
ing \Ge\ = IGa/I, the predicted branching ratios are 
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FIG. 4. The differential decay rate vs t = rrip^ of the vec- 
tor current induced -B" — > D*~pn decay. The upper curves 
are from the phenomenological fits to nucleon form factor 
data assuming \Ge\ = |Ga/|, for the Light-Front (solid) and 
BSW (dashed) B" -^ D*^ form factor models. The lower 
curves are for the VMD model discussed in the next section. 



are smaller than the BSW model case. The ~ 10% dif- 
ference can be viewed as an estimate of the uncertainty 
from B -^ D* form factor models. 

From Fig. Q we sec that the differential rate peaks at 
~ 4.6 GeV^, corresponding to riipn — 2.14 GeV, which is 
quite close to the threshold of 1.88 GeV. This threshold 
enhancement effect, consistent with what was suggested 
in Ref . ^ , should be checked experimentally by measur- 
ing the recoil D*^ momentum spectrum. 



D. Estimate of Axial Current Contribution 

Although the time-like data for the form factors of the 
axial current is still lacking, we can nevertheless make a 
rough estimate of its contribution. In analogy with the 
nucleon em form factors which are constrained by the 
asymptotic form of Eq. (po|), we expect that, for large i, 
gA (t) behaves as 1/i^ and dominates over Ha (i), which 
behaves like 1 /t'^ . Taking cue from the similarity be- 
tween Eqs. (|^ and (|lj), we estimate the axial vector 
contribution by making a simple comparison and scaling 
from the vector case. Since we only have space-like in- 
formation for the gA{t) form factor, we estimate gA{tth) 
at threshold by assuming same threshold enhancement 
factor as the G^f{t) — C^.At) case. 

With this and Eqs. (17) and (18) in mind, to es- 
timate the axial vector current contribution, we scale 



the decay rate obtained from the vector case by the ra- 
tio gi(-ith)/(G^/(-ith)-G^(-ith))^ for space-like mo- 
menum. We use a dipole fit gA{t) — 5^(0)/(l — t/M'j)'^ 
with the axial mass Ma = 1.077 ± 0.039 GeV H. For 
t = -4m^, the ratio r(i) = gA{t) /{GP{t)-Glj{t)) gives 
r(— 4™^) ~ 0.59. Assuming this ratio holds also at the 
threshold t — Amj^, then r^(4m^) ~ 0.35 could be taken 
as the ratio of the branching fraction from the weak axial 
vector current to that from the weak vector current, i.e. 
Bi a/Btv ~ r^(4m^). For the results from fitting data 
assuming \Ge\ — \Gm\, the total rate Br = Bry + Br a 
would then be Br^^ ^ {I + 0.35) x (7.14 x lO"") = 
9.64 X lO"'' and Br^^^ ~ 11.77 x lO""*. On the other 
hand, for the results from fitting data assuming \Ge\ = 0, 
the total rate would be Br^'' - (l-h0.35)x (8.96x10"'') = 
12.10 X 10-4 and Br^^^ - 14.77 x 10"^, quite compatible 
with the experimental value of (14.5 ^gg ± 2.7) x 10"''. 
Another value of t which could be of interest is i = 0. 
Assuming that the ratio r(0) ^ 0.27 holds at thresh- 
old, then from r^{0) ^ 0.07 ^ Bta/Btv, the total rate 
Br = Bry + Bta is dominated by the vector current con- 
tribution. 

To improve our result, we urge for more experimental 
measurements of G^j(t), G^^"'{t) and gA{t)- In turn, if 
the predictions (strength and spectrum) from our model 
based on the vector part are confirmed by experiment, 
the meaurements of B — > D*pn could provide useful feed- 
back on nucleon axial form factors. 



III. INTERRELATION WITH NUCLEON FORM 
FACTOR MODELS 

The nucleon form factor is one of the oldest subjects 
in particle physics. It has provided us with a wealth 
of information and insight, and remains an active field 
to this date. In the preceding section, we made a simple 
phenomenological fit of nucleon em form factor data, then 
used an isospin relation and a factorization hypothesis 
to compute the i?° -^ D*^pn rate, with some success. 
Although we made use of PQCD counting rules, we did 
not utilize tools such as analyticity. On the other hand, 
we mentioned the possibility that B'^ — > D*~pn type of 
modes could eventually provide useful information on the 
nucleon form factor itself. 

To exploit the utility of analyticity and to illustrate 
future interrelations between B'^ — > D*^pn and nucleon 
form factors, we adopt a specific nucleon form factor 
model and discuss where it may be improved. The dis- 
cussion would also shed some light on form factor decom- 
positions, as well as the possible approach to i? — > 7pA, 
■q'pK modes, which we shall briefly touch upon later. 



A. Dispersion Analysis and VMD Hypothesis 

Among the various approaches to the nucleon em 
form factors, of particular interest is the Vector-Meson- 
Dominance (VMD) hypothesis, which states that a pho- 
ton couples to the hadrons via intermediate vector 
mesons such as p, w, 0, etc. The simple pole model 
of Fig. |l|(a ) is a limited form of the VMD hypothesis. 
In Ref. p2| , a parameterization of the nucleon em form 
factors based on dispersion analysis was proposed, which 
is constrained by several physical conditions, including 
PQCD power counting. The starting point is the disper- 
sion relation 



F{t) = - 

IT 



to 



Im F{t') 
t' -t 



dt' 



(31) 



where F{t) stands for the nucleon em form factors F^2{t)- 
To gain predictive power, one truncates the spectral func- 
tion by a few vector meson poles, where, to mimic the 
effect of large t continuum and to enforce PQCD count- 
ing rules, a fictitious pole is introduced. Thus, the form 
factors take the form 



F:{t)^Y. 



M?-V 



1,2, 



(32) 



where a\ ^ are related to the vector (« = 1) and ten- 
sor (i = 2) coupling constants of Eq. (|^) via 



V2< = m^!vgl 



vNN 



(33) 



and /i, is the vector meson decay constant. This clearly 
extends the simple p meson exchange picture. Together 
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TABLE I. The relevant poles (Af„ in GeV) and the corre- 
sponding residues that enter the summation in Eq. (b3) 



M,,^ 


1.4035 


M,„ 


= 1.45 


Mp,„ -- 


= 1.69 


< 


"2 


< 


0-2 


< 


p'" 

0-2 


-9.913 
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FIG. 5. The VMD picture for B° --> D*-pn decay. The W 
boson couples to the nucleon pair via vector mesons which in 
our case are p and three excited states, including a fictitious 
pole to mimic continuum efi'ect. 



with our factorization ansatz, the i?*^ —^ D*~'pn transi- 
tion can be pictured as in Fig. g. 

Following Ref. |2^ we now consider the following pa- 
rameterization of the form factors: 



Fr(t) 



^"w+E 



<L-\Ml) 



M^-t 



Lit), 



where 



m^ 



In 






(34) 



(35) 



Ftit) 



a'iL-\Ml)+b<lL-\Ml){l-t/c<[) 
2(1 -V<) 



-2/i 



(36) 



4 



for i = 1, 2, where a^ = 1.0317, a^ = 5.7824, b^ = 0.0875, 
0.3907, c^* = 0.3176, 4 = 0.1422, < = 0.5496, 
0.5362, and M^ = 0.5 GeV^, M^ = 0.4 GeV^. 
Since the form factors also receive constraints from per- 
turbative QCD for large momentum transfer, a logarith- 
mic factor, Eq. (p5|), is given for consistency. Apart from 
this, Eq. ( |34| ) contains two terms: F[{t) represents the 
27r-continuum plus the p pole, the remainder a summa- 
tion over additional isovector vector meson poles. 

The VMD model of ref. p^] focused more on fitting 
the space-like nucleon form factors. It was found suf- 
ficient to use 3 additional vector meson poles, two of 
which are chosen to be the physical particles p(1450) and 
p(1700) and denoted as p" and p'" . In Ref. [l|], which 
extends the scenario to include time-like data, the third 
pole is left adjustable to compensate for the neglect of 
higher mass continuum like NN. The pole mass was de- 
termined to be Mp' = 1.4035 GeV, in association with 
fixing A2 = 12 GeV^ and Q^ = 0.35 GeV^ in Eq. ^. 
The parameters in Eq. ( p6|) remain unaffected in both 
fits. Besides the fictitious p' , one special feature of the 
model is to utilize the freedom in insufficient knowledge 
of p" and p'" widths and couplings, which are taken as 



summation so that Ff{t) contains no parameters that 
need to be determined. 

With a^ taken as parameters, Eq. (|2C|), i.e. PQCD 
power counting, can be achieved by choosing the residues 
a^ of the vector meson poles such that the leading coef- 
ficients in the 1/i expansion cancel. This is what the use 
of a single p-pole can never achieve, since further excited 
states are needed for such cancellation. This also means 
that one has only an effective model since the dynami- 
cal parameters for p{lA50) and p(1700) would likely not 
correspond to physical ones. We summarize in Table m 
the relevant meson poles and the corresponding residues 
of the higher excited states given by Table 1 of Ref. |l^] 
(only the so-called "Fit 2" is needed). 

The parametrization agrees with experiment quite well 
up to i ^ 6 GeV^ , beyond which it overruns data because 
of the choice of low A in the logarithmic factor L{t). This 
is in contrast with the empirical fit of Eq. (|2l|). It arises 
in part because the VMD model focuses more on the 
spacelike region where one has more data, but is of less 
concern to us. In order to conform to experiment for 
larger timelike momentum transfer, however, a modifi- 
cation of the L{t) factor is needed. The empirical fit of 
Eq. (pi]) suggests a convenient modification 



Lit) 



In 



A2 
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Qi 





1 Tl ( A'' 
"2A^ \~ 
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Hm 



A2 

for t < A, 

- 2 

^y+H{A) 

for |i-^|<A, 

A2 
for t> h A. 



where we match a parabola between the interval ^ — A < 
t < ^ -t- A by tuning the constant H{A). Note that 
-L' (AV2 - A) = L' (AV2 + A) and L {K^ /2 - A) = 
L [h? /2 + A), or L{t) is symmetric with respect to t — 
A^/2. To smooth out the artificial rise that occurs be- 
yond t ~ 6 GeV^ = A^/2 for the original fit, we choose 
A = 0.5 GeV^ and H (A) ~ 0.10279. 

Fig. II shows the resulting VMD-based proton mag- 
netic form factor for both space-like and time-like i, with 
the modified log factor in the time-like region, where the 
result is plotted in more detail in Fig. ||. We see that 
the trend of the proton data can be described by this 
model. In contrast, from Fig. 0, where the extension of 
the VMD result to neutron case is given, there is a signif- 
fit parameters. Thus only the p(770) is isolated from the icant deviation from G^ data, signaling the incapability 
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FIG. 6. The VMD-based proton magnetic form factor 
with the modified logarithmic factor. The unphysical region 
is marked by two vertical lines. The dash-line in the time-like 
region representsa fit with C = 53 in Eq. (|2ll). 



of the thrcc-plus-one pole fit to describe the neutron data 
consistently. Thiis fact was admitted to in Ref. p^ . 

One useful aspect of a dispersion approach is that ana- 
lyticity can help determine the signs of the time-like form 
factors from space-like region. Unlike Eqs. ( p^ ) and ( |2^ ) 
where the signs are put in by hand, it is more natural in 
the dispersion analysis since all the time-like information 
can in principle be obtained via the dispersion relation 
in Eq. (pTj), where the spectral function Iva F{t) is ana- 
lytically continued from the space-like region. One can 
readily check this by finding out the value of the mag- 
netic form factors from the VMD analysis at threshold: 
Gltiimjf) ^ -fO.39 while Cljiimjf) '^ -0.22. From 
both Figs. H and & 
cross the i-axis, G^j 



since neither G^^ 



nor G^ seem to 



remains positive while G^.j- remains 
negative throughout the timelike region. VMD analysis 
thus gives G^j and G^j with a relative sign. 



B. B° ^ D*-pn in VMD Approach 
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FIG. 7. C^j — G'm (long-dashed) from our fit to nucleon 
form factor data compared with 2Fi + 2F2 {— G\j — GJ/ , dot- 
ted) from the VMD analysis, for the kinematically allowed re- 
gion oft in B^ — > D*~pn. Also shown are 2Fi (short-dashed) 
and 2F2" (dot-dashed) from the VMD model. 



for both LF and BSW models. Had we used Eq. (|3£ 
because of the artificial rise in this original logarithmic 
factor, the contribution from t > 6 GeV^ {mpn > 2.45 
GeV) would be ~ 2.5 times higher. The contribution 
from t < 6 GeV^, however, is unchanged. 

It is clear that the branching fractions obtained in the 
VMD approach of Ref. m& are typically 5 times smaller 
than our phenomenological model discussed in the previ- 
ous section. This can be simply traced to the inadequacy 
in accounting for neutron data by the VMD model, as 
seen by comparing Figs. and pi While giving a rea- 
sonable fit in the proton case, the absolute value of its 
neutron form factor is simply below all data points. The 
VMD model was tuned more on the proton where data 
is much more abundant. Since the proton and neutron 
magnetic form factors are opposite in sign, if the VMD 
approach is improved to give better account of |G"^ | data, 
the combination Gl'j — G^j would increase. 



We calculate the branching fraction of B^ -^ D*^pn 
through the vector portion of the charged weak current 
in the VMD model, again taking ai as extracted from 
B" -^ D*-p+. The results for S" -^ D* 



pn are. 



BrV/ 



(1.69 X 10""*) ( 



0.74 



Br 



BSW 
V 



(2.06 X 10" 



^0.86 



(37) 
(38) 



which are respectively about 12% and 14% of the exper- 
imental value of (14.5 +i^o ± 2.7) x 10""^. Varying pole 
masses slightly does not drastically modify the result. 
The differential decay rate is plotted in Fig. H, where 
one sees that it peaks not far above the pn threshold 
oi t ^ 3.53 GeV^. The contribution from the region of 
3.53 < t < 6 GeV^ is more than 80% of the total rate 



IV. DISCUSSION AND APPLICATIONS 



A. Nucleon Form Factor Decompositions 



We plot various form factor combinations in Fig. |^. 
The long dashed line shows GIj-GIj = 2{F^ + F:^) from 
our phenomenological fit, which assumes G^j and G^li 
have opposite sign. As discussed in the previous section, 
because the VMD model gives much lower value for | G^j \ , 
G^,^ — G^ in the VMD model (denoted as dotted line and 
labeled by 2{Fi + F2)) stays below the phenomenological 
model. However, had we chosen the proton and neutron 
form factors to be of the same sign , G^^ — G^^j for our 
phenomenological fit would be very close to the solid line 
which is the t-axis, and would give a rate that is two 
orders of magnitude too small. 



TABLE II. Comparison of nucleon form factor decomposi- 
tions. The fraction TZx = Sx/Brv are defined in text, where 
X stands for any combination of _Fi, F2 or their product. 





^(F,"+F,") 


TZ,. 


^(F,''+F,")-F,'' 


LF/BSW 


1.3 


0.24 


-0.54 




7^F, 


7^F, 


7^F^F," 


LF/BSW 


26.00 


18.25 


-43.25 



Besides helping to fix the sign of G^^ by analytic- 
ity, another utiUty of discussing the VMD approach is 
that it can give some insight to the F2 nucleon form 
factor. Because of lack of experimental information, we 
have dropped the F2 contribution in our phenomenolog- 
ical approach, and we need to check the validity of this. 
The weak vector current induced decay amplitude, upon 
squaring, can be expressed as 



\Mvf =\MF,+Ff 



\MF,r 

+2Re {Mf,+f,M*fJ 



where MF1+F2.F2 denote F^ + F2 , F^ terms in Eq. (|17[). 
Decomposing Bry = Bf^+f^ + ^f^ + ^(Fi+F2)F2 where 
the last term is the interference term, we define the rel- 
ative fractions such as TZ(^Fi+F2) = ^(Fi+F2)/^'^v from 
{F^+F^) alone. We find 7^(f,+f,),7^F,,7^(f,+J.,).f, = 
130%, 24%, -54% (Table 11) in the VMD model, for both 
LF and BSW B° -^ D* form factor models. Note that 



the 2(F]" -f F^ 






Gn 
M 



term gives the dominant 



contribution, which supports the approximation used in 
Sec. II. The Fj contribution is much smaller even though 
1^2^ (i) I is greater than \F^ + F^\ for t < 5 GeV^, as shown 
by the dot-dash curve in Fig. [7[ The interference term 
contributes ~ 40% of the F^ + F^ contribution. The 
destructive nature is due to the relative sign between 
Fi + F2 and i^^, which reduces the combined effect of 
the last two terms in Eq. (|39|). 

It is instructive to compare with Eq. (0), where one 
decomposes into Fi and F2 directly. As shown in Ta- 
ble O, the individual terms from this decomposition are 
an order of magnitude larger, and only after strong can- 
cellations does one arrive at Bry. It is therefore not a 
very useful decomposition. We see from Fig. R that the 
magnitudes of F^ and F2 are all greater than their sum, 
hence F^ + F^ together with F^ is a better decomposi- 
tion for t close to threshold, which is what we used in our 
phenomenological study. 



B. Prediction for B -^ D^*^pn Modes 

Our phenomenological approach can be applied to the 
modes of i3+ — > D^*^"ph and B^ -^ D~pn. We show 
in Table |l| the results based on the vector current with 
the B'^ — > D*^pn mode included for comparison. The 



TABLE III. Branching fractions of the B -^ D^*'pn modes 
from the vector current, obtained by using the phenomeno- 
logical form factors with \Gl.j\ = |G^| and \Gli\ = \Ge\ 
{\Ge\ = 0) for the first (second) set. 



Brv X 10* 

\G"m\ = \G%\ 
LF BSW 



Brv X 10" 
|GSI=0 
LF BSW 



B+ 



D* 
> D' 



pn 
"pn 



7.14 
7.64 



B+ 
B° 



+U.69 

0.65 

+0.74 

-0.69 



8.72 

9.33 



> D^pn 
D^pn 



3.92 
3.66 



+u.3y 
0.35 

+0.36 
0.32 



+U.85 

0.79 

+0.90 

-0.85 



8.96 
9.59 



+ LW2 



1+1.^5 



0.94 l'-'.94_j^ J5 



+ 1.10 



-1.33 



3.21 
2.99 



+U.3Li 
0.28 

+0.30 
0.26 



4.89t-- 4.00t^- 

-+O.54 o yQ+0.44 

-0.48 ''•''J-0.39 



4.56! 



differential decay rates for the other three B -^ D'^*^pn 
decays are given in Fig. g. We have used the central val- 
ues of the effective coefficients that are extracted from 



the two-body modes [E§ 5° ^ D 



*+, 



D'^ p with val- 



ues a 
and a 



LF(BSW) 



LF{BSW) 



0.74 ± 0.18 ± 0.06 (0.86 ± 0.21 ± 0.07) 

0.89 ± 0.08 ±0.07 (0.91 ± 0.08 ±0.07) 

are 



(39) respectively. We note that the D* modes in Table [II 



close in rate, and likewise for D modes, which is easy to 
understand because of simple replacement of spectator 
quark in B — > D^*' transition. 

Although only the contribution from the weak vector 
current can so far be calculated, we can estimate the 
values of these other baryonic modes by following the 
recipe of Sec. II. Inspection of Table III shows that the 
ratio of Brv/(B+ -^ D*°pnyBTv{B° -^ D*-pn) remains 
fixed regardless of S ^ D* form factor models. As- 
suming similar behavior for axial contribution, we expect 
the B^ -^ D*'^pn branching fraction to be scaled by the 
same factor and find the value of ^ 15.5 x 10^'* that 
is only slightly larger than the i3° -^ D*^pn mode, as 
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FIG. 8. The differential decay rates arising from nucleon 
vector current. Upper (lower) dashed line is for B" — » D'^^pn 
with the BSW (LF) hadronic form factors; upper (lower) dot- 
ted and solid lines are for the B" -^ D'^pn and B" — > D^pn 
modes using the LF (BSW) model. 



TABLE IV. 
BiviB -^ D^*^pn)/'B,vv{B° -^ D'-pn) of Table |I| with re- 
spect to the Bt{B" -^ D*~pn) experimental central value of 
14.5 X 10-'' (Eq. (1)). 



1 bv 





BrxlO* 


BrxW 




\r<" 1 — ir^" 1 
l"j"A/l — I^-^bI 


\G"e\=0 




LF BSW 


LF BSW 


B+ -^ D*"pn 


15.52 15.51 


15.52 15.51 


B+ -> D^'pn 


7.97 5.34 


7.92 5.30 


B° -^ D~pn 


7.43 4.98 


7.38 4.94 



given in Table |rV|. The predicted branching fractions for 
5+ — > D'^pn and B'^ -^ D~pn from the same ansatz are 
in general are 2-3 times smaller. Since one is comparing 
B ^ D vs. B ^ D* , there is stronger model dependence 
on the transition form factor. 



One can see from both Tables II] and IV as well as 
Fig. ^ that the LF results are generally larger than the 
BSW ones ior B ^- D modes, but the case is reversed 
for the B -^ D* modes. This can be understood from 
the differences in ai and hadronic form factors of the 
LF and BSW cases. For the B ^ D modes, the only 
hadronic form factor involved is F^^, which behaves as 
a dipole and monopole in the LF and BSW models, re- 
spectively. This leads to F^-^it) > Ff^'^{t) in the phys- 



= 0.9, giving a 



ical timelike region while af^ = af^'^ 
larger Br^^. In contrast, for B -^ D* modes, although 
the dominant hadronic form factor A^^ > Af^^ in the 
physical timelike region, it behaves as monopole for both 
LF and BSW models, hence the difference between Af ^ 
and Af^'^ is not as large as the previous case. However, 
ai^ ~ 0.74 < af^^ ^ 0.86 hence we have the opposite 
result that BSW rates are bigger. 
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FIG. 9. Ff + F^ for the pn (dotted), E°E- (solid), 
H^H" (dashed) and S^A (dot-dashed) modes, respectively. 
The vertical lines specify the thresholds for the baryonic pairs. 



We shall only give results from the phenomenological ap- 
proach with \Gl\ = \Gli\. 

The strange baryon modes all have rates smaller than 
that of the pn mode due to the following reason. In 



w 



Fig. 1^ we plot the form factor combination F^ + F. 
2 (F^" -1-^2") of Eq. (|l^. One can see that the largest value 
is at the pn threshold, while for the strange baryonic 
modes the corresponding threshold values are all smaller. 
Reading off from Fig. 0, it is clear that the E+A mode 
would be the dominant strange baryonic mode, with the 
two SS modes the smallest. This is shown in Fig. HG 
for the differential decay rates using LF mesonic form 
factors, and the total decay rates given in Table M. Note 
that the differential decay rate for the Z)*-Eli mode has 
a zero at t ~ 7.5 GeV^ because F^(j:t) + F^{J:t) 
changes sign, as can be seen from Fig. \ 



C. Predictions for B'^ — > D* + Strange Baryons 



Our phenomenological approach can be extended to B^ 
decay into D*^ plus baryon pairs containing strangeness. 
Recall that in Sec. O we utilized SU(2) symmetry to ob- 
tain the relation F{^ = -PT2 ~ -^i"2 fo'" the pn mode. In 
the SU(3) limit we can |24| extend this relation to modes 
containing strange baryons. Starting from Eq. (^, de- 
noting F]'^(BsB^) as the weak form factor that appears 
in the matrix element (BsBj|V"+|0), we find 

FW(^^+,oj.o,^) = tV2 (Fl.it) + i F-,{t)] , 



F^,iE"E-)=Fl,it) + 2Fl[,{t), 
F^,{i:+A)^-J^Fl],{t). 



(40) 



These relations enable us to calculate the contributions 
from the vector currents to the strange baryonic modes. 
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FIG. 10. Differential 

rates of B° -^ L>*-E+A (dot-dashed), 1)*-=°=- (dashed) 
and D*~'S^TP (solid) from the phenomenological model, with 
LF hadronic form factors. 
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TABLE V. Branching fractions for B° -^ D*^ + strange 
baryon pairs from vector current using the phenomenological 
Gm form factors of Sec. Rassuming |GJ| ~ \G^j\. 





B 


r\r X 10* 


Br^sw X 10* 


B" -^ D*-pn 




7.19 


8.78 


B° ^D*-E+A 




1.09 


1.38 


B" -^ D*-H°E- 




0.52 


0.66 


B" -^ D*-E+E° 




0.01 


0.02 


B" -^D*-E°E- 




0.01 


0.02 



V. CONCLUSION 

The main result of this paper is an attempt to ac- 
count for CLEO's result on B^ -^ D*~pn. In a fac- 
torization approach, the nucleon pair is viewed as pro- 
duced directly from the weak current. We then use an 
isospin relation of weak and em vector currents to ob- 
tain the vector weak nucleon form factors directly from 
their em partners, where a relatively large database ex- 
ists. It is interesting that we can account for ~ 60% of 
the observed Br(i3" -^ D*~pn) rate in this way. A VMD 
model that attempts at fitting nucleon form factor data 
was discussed to clarify certain issues. 

Interference of the weak vector current induced am- 
plitude with the amplitude induced by the weak axial 
current does not manifest itself in the total rate. The 
latter is a simple sum of the absolute squares of both 
vector and axial vector. A rough estimate of axial vector 
contribution, together with the dominant vector contri- 
bution, seem to fit the measured rate. However, for a 
more reliable calculation, more measurements on time- 
like region G^" and gA are urged. 

We apply our analysis to B+ -^ I)^*^^pn and i?° -^ 
D~pn modes to get the rates arising from the vector 
current. We find Br{B+ -> D*^pn) ~ Br(B° -^ D*-pn) 
and Br(B+ -^ D°pn) - By{B° -^ D~pn), with the lat- 
ter modes having smaller rates slightly below the 10"'^ 
level. Our analysis is also applied to baryonic modes 
that contain strangeness. The estimated branching frac- 
tions are generally lower than the pn mode due to smaller 
couplings and higher thresholds. The largest modes, 
B" -^ £)*-!]+ A, is predicted at the 1 x 10"'* level. 

For the analogous picture of i? -^ v'^P^ l^P as de- 
scending from B -^ rj'K and K*j via K^*^ exchange, 
the baryon pairs are not produced by the W boson, and 
it seems that the approach used here cannot be read- 
ily applied. However, some experience obtained may 
still be valuable. For example, in the VMD approach 
to i? — > D*pn, more than one pole and cancellations 
among them are required to reproduce the QCD pre- 
dicted asymptotic behavior. For the charmless cases, the 
baryon pairs are no longer produced by the charged cur- 
rent. Instead, the K'^*'> resonances that correspond to 
pi,ii.Jii jjj VMD approach, appear more as string excita- 



tions. They need not obey the same QCD power count- 
ings, and perhaps may result in larger rates. 

Finally, it is of great interest to estimate the rate of 
the charmless baryonic mode B" — > p~pn by replacing 
D*^ in the Feynman diagram depicted in Fig. 0(b) with 
p~ . Since this is a tree-dominant mode, extending from 
the study presented in this paper, Br(i?*^ —^ p~pn) could 
be as large as that of the two-body mode B^ -^ P~ P^ , 
analogous to the relative strength of Br(i?° -^ D*~pn) 
vs Br(i3° ^ D*~p+). Estimating via Br{B° -^ p~pn) ^ 
Bt{B° -^ p^p+)x Bt{B° -^ D*-pn)/BT{B° -^ D*' p+). 
From Bt{B° ^ p- p+) - (20-40) xlO'^ 0,|| we get 
Bt{B° -^ p-pn) - [0.4-0.8) xlO'^ Alternatively, we 
can scale from B -^ D*pn by \Vub/ycd\^ and phase space, 
decay constant etc., and again find B — > p^ph at 10~^ 
order. Charmless decays are of great current interest. 
A more detailed discussion of i? ^ ppn modes is given 
elsewhere lEQl . 
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